Coexistence of ferromagnetism and Kondo effect in uranium compounds 
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Coexistence between ferromagnetic order and Kondo behavior has been observed in some uranium com- 
pounds. The underscreened Kondo lattice model can provide a possible description of this coexistence. Here we 
present a model of a lattice of 5 = 1 spins coupled to the conduction electrons through an intra-site exchange 
interaction Jk and an inter-site ferromagnetic exchange f — f interaction Jh ■ Finite temperature results show 
that the Kondo temperature is larger than the Curie ordering temperature, Tc, providing a possible scenario 
for the coexistence of Kondo effect and magnetic order. Also, the Kondo behavior disappears abruptly for low 
values of Jk and smoothly when changing the band occupation. These results are in qualitative agreement with 
the experimental situation for the above mentioned uranium compounds. 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr, 75.10.-b 



Many rare-earth and actinide compounds show anomalous 
properties which have been attributed to heavy fermion behav- 
ior. A lot of attention has been placed on the experimental and 
theoretical study of the competition between magnetic order - 
mainly antiferromagnetic (AR - and Kondo effect in cerium 
or ytterbium compounds As a result, nowadays there 

exist various available theories describing many reliable ex- 
periments, or vice versa. For a review see ref. Li]- 

On the other hand, uranium compounds present many dis- 
tintic behaviors from mixed valence to magnetic order with or 
without a Kondo effect and even a coexistence between feiTO- 
magnetic order and superconductivity, as it has been observed 
in UGe2 and URhGe ll. 

However, one of the most interesting experimental find- 
ings, the coexistence of Kondo behavior and feiTomagnetic 
order observed in some uranium compounds has been - to our 
knowledge - somewhat overlooked from a theoretical point of 
view. It's a long way since the first experimental evidence of 
coexistence between Kondo behavior and ferromagnetic or- 
der in the dense Kondo compound C/Te yj. More recently 
this coexistence has been observed in UCuo.gSb2 |5] and 
UCoQ,5Sb2 |6]. Those systems undergo a feiTomagnetic or- 
dering at relatively high Curie temperatures of Tc = 102K, Tc 
= 11 3K and Tc = 64. 5K, coiTespondingly, and also present 
a logarithmic decrease of the magnetic resistivity above it, 
suggesting a Kondo behavior at high temperature. These ex- 
perimental results clearly show the coexistence of feiTomag- 
netic order and Kondo behavior below a large Curie temper- 
ature. Another example of such coexistence has been found 
in U Ru2-xRe.xSi2 compounds, where a non-Fermi liquid 
(NFL) behavior of resistivity and specific heat has been ob- 
served inside the ferromagnetic phase |7, 8]. 

In order to account for the Kondo-ferromagnetism coexis- 
tence observed in the above quoted uranium compounds, we 
propose the framework of the underscreened Kondo lattice 
(UKL) model, that incorporates all the essential aspects of the 
problem such as Kondo interaction and magnetic RKKY in- 
teraction between localized spins = 1. We believe it is 



reasonable to assume a bp configuration for uranium because 
the magnetic moments deduced from magnetic experiments in 
these compounds are close to the free ion values 1 9] and more- 
over a good agreement with experiments has been obtained by 
considering only 2 localized 5f electrons 1 10]. 

The underscreened Kondo impurity model has been studied 
using different approaches fl7| and an exact solution was ob- 
tained by the Bethe-ansatz method |12]. On the other hand 
there are relatively few studies of the UKL, among them 
refs.|13, 14]. Here we describe the UKL as a periodic lat- 
tice of magnetic atoms with S = 1 - which can be modeled 
by two degenerate /-orbitals - interacting with itinerant spins. 
The spins of these /-electrons are coupled to S* = 1 due to the 
strong on-site Hund's coupling and we treat them on the ba- 
sis of a fermionic representation for spin S = 1 ]15], valid 
in this constrained triplet Hilbert space. Then, we perform 
a mean field analysis of the model, both at zero- and finite- 
temperatures in terms of the Green's functions of the system. 

In normal Kondo lattice (KL) model {S = 1/2, i.e. one 
/-electron per lattice site) the competition between Kondo ef- 
fect and magnetic RKKY interaction leads to the appearance 
of a quantum phase transition between a magnetically ordered 
(generally antiferromagnetic) phase and a phase with coher- 
ent Kondo spin-singlet formation and short range magnetic 
correlations In a similar way, with the help of the 

UKL model we can describe quantum feiTomagnetic phase 
transitions and the coexistence between ferromagnetic order 
and NFL behavior ]|8j] at different values of external parame- 
ters such as band filling, pressure or temperature. Based on 
this analysis we can draw a phase diagram showing magnetic 
ordered phases with no Kondo effect and regions of coexis- 
tence between ferromagnetic order and Kondo behavior 

The UKL Hamiltonian can be written in the following form: 
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Jk ^ SiCTi + -Jh ^ SiSj 



(1) 



where the first term represents the conduction band with dis- 
persion €k, width 2D and a constant density of states 1/2D, 
while /i is the bare electron chemical potential. Localized 
spins are represented by fermionic operators f^^^ and iiaa 
introduced in refs.|12, J^, carrying spin, a, and orbital, a, 
indexes. Eq is a Lagrange multiplier which is fixed by a 
constraint for the total number of f-electrons per site, n/ = 
X](t("-o-^ + n^^) = 2, and can be interpreted as a fictitious 
chemical potential for /-fermions. The third term is the on- 
site Kondo coupling, Jk > 0, between localized Si = 1 and 
conduction electron's ai = 1/2 spins, whose fermionic repre- 
sentation is done in the usual way. And finally, the last term 
is the ferromagnetic inter-site interaction, Jh < , between 
localized /-magnetic moments, resulting explicitly from two 
contributions: the effective RKKY interaction, and the direct 
exchange. 

As a first step we introduce the relevant fields. They are 
AjCT = I]„ f« which couple Ci^ and electrons at the 
same site, and operators of magnetization for both c- and /- 
subsystems, which are, respectively. Mi — Sf = \{nl^—n^^) 
and riii = af = ^ (rt^^ — ) describing the long range mag- 
netic order of the system. Then we perform a mean field (MF) 
analysis of the Hamiltonian eq. (1), expanding it according 
to the four following averages: = {Xia), M = (Mi) and 
771 = (mi). The non-zero values of M and ttj correspond 
to the formation of the magnetic phase with finite total mag- 
netization, while the parameter describes an effective hy- 
bridization between conduction and /— electrons that corre- 
sponds to the Kondo behavior In MF approximation the aver- 
age value of c^afia be different from zero, implying 
spurious charge fluctuations. But one should keep in mind that 
the original term of the Hamiltonian is a four-fermion opera- 
tor, i.e. a product Ao-A^ 1 1?]- So, Ac, should be interpreted as 
a tool for the MF calculation that provides a correct descrip- 
tion of the Kondo temperature, for example in the study of the 
Kondo impurity |18]. This procedure is basically equivalent 
to other mean field approaches developed for studying nor- 
mal KL, such as the path-integral calculation restricted to a 
saddle-point solution performed by Coleman and Andrei 1 1^, 
or large-A^ formulation performed by Burdin et al |20], where 
again the saddle-point solution yields a mean-boson-field ap- 
proximation. 

The resulting MF Hamiltonian reads 



HmF - J2ka ^ka'^la + Siaa EoalT-taa^ 



(2) 



2JKNJ2a AffA^ - JkNtuM - \JhNzIVP 
where, remembering that Jk > and Jh < 0, energies read 



e^^ — £j: + Jk<^M , with a = ±1/2 
Eoa = Ef) + JKcrm - ^JK^a^a + JhzctM 



(3) 



The diagonalization of the MF Hamiltonian (|3 gives two 
non-hybridized /-bands (one for each spin) at energies Eqo- 
and two quasiparticle bands Ej.{k) with energies 



Em = -[Eo. + sj:^ ± ^{Eo. - Sj^^f + 8a|] (4) 

where ag ~ ~\JkXs- The ± sign refers to the upper (lower) 
hybridized band. Let us point out here one of the main differ- 
ences between the mean field treatment of the normal KL and 
UKL models: in the underscreened case one /-localized level 
remains non-hybridized and for the other level the resulting 
c — / effective hybridization is twice that of the normal case. 
The energy spectra E'^ (k) depend on a set of external parame- 
ters such as band filling Uc, Kondo coupling Jk and exchange 
interaction Jh and a set of internal parameters, M, m, X^, fi 
and Eq, which should be calculated in a self-consistent way. 
The system of self-consistent equations can be obtained by 
keeping constant the number of /- and c- electrons. Their 
expression can be evaluated from the Green functions of the 
system by straightforward calculations : 



'V ~ 2D J- 

^f{E+„Y-^^-^ 



de„ [nF{Eo„)- 
+ nAE^.Y-^] 



"c 2D J- 

npiE 



D+A^ 

2D J-D+A„ 



d£c,[~nF{E+a) 



En. 



(5) 



where np{uj) 



— -typ 



is the Fermi distribution function. 



= JkctM and Wa{e) = y/iEoa - e)^ + 8a|. From 
the Green functions we can also obtain the expression for Acr 
given by: 



- h !- 



D+A„ 



D J-D+A, 



dea[nF{E+a) - 71f(£;_^)]t^ 



(6) 



In order to construct a set of self-consistent equations de- 
fined by (|5} and (|6j, we substitute the strong constraint, riij — 
2 in each site, by a softer one for its average, and set the band 
filling by the averaged number of c-electrons, ric, which we 
consider in the usual range of partial filling, i.e. < 77c < 1 : 



Uf = 77 I + 7li ~2 



(7) 



Having solved the system of self-consistent equations and 
also minimized the free energy, we study various properties 
of the UKL model. In Fig Q we present a summary of the re- 
sults for r = 0. The region of the parameters Jk and 77c 
that provide finite values of A| corresponds to a coexistence 
between the heavy fermion behavior and ferromagnetic order 
The values of A| are close to A| so the former is not depicted 
on the figure, neither the magnetization, to preserve clearness. 
It is possible to see on Fig^that A| decreases smoothly as a 
function of ric, following an approximate square root behav- 
ior J3l, while it undergoes a sharp transition as a function of 
Jk- When A| goes to zero the ground state is magnetically or- 
dered with no Kondo effect. The hybridization gap — 8a? 
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FIG. 1: Plot of the parameter A| as a function of Jk and ric for 
T = 0. The ferromagnetic interaction between localized /-moments 
is Jh ~ —0.01. There is a discontinuous transition as a function of 
Jk and a behavior in y/rh- as explained in the text 



vanish simultaneously, signaling a quantum phase transition. 
Detailed calculations and results including a full phase dia- 
gram with different values of Jk and Jh will be presented in 
a forthcoming article. 

Finite temperature calculation permits to determine the 
Curie temperature and the correlation (Kondo) temperature. 
Fig|2l shows the behavior of the /— and c— magnetizations, 
M and m respectively, and A^, as a function of the tem- 
perature. The parameters here and in following figures are 
Jk — 0.8, Jh = —0.01 and ric = 0.8 and all energies and 
temperatures are in units of the half -bandwidth D. At zero and 
low temperatures we observe the coexistence of magnetic or- 
der and heavy fermion behavior, but as both spin systems are 
strongly polarized the Kondo effect is concealed. When the 
magnetization decreases, A^. grows, having its maximum at 
the Curie temperature when the magnetization vanishes. Due 
to the breakdown of the spin symmetry A| and A| are slightly 
different in the magnetic region but they coincide when the 
magnetization vanishes for T ~ Tc, i.e. when the spin sym- 
metry is restored. For T > Tc the system exhibits only Kondo 
behavior (Acr 0, M — and m ~ 0). Finally, above a char- 
acteristic temperature, the Kondo temperature, Tk, both A^ 
vanish and the two electron systems are decoupled. 

This behavior can be further clarified by the plot of the den- 
sities of states at different temperatures. As the c— densities of 
states are almost constant (except in the hybridization gap re- 
gion) we present the results only for the /—densities of states, 
which are calculated numerically from the imaginary part of 
the / — / Green functions. In Fig.Q we plot /0/cr(e — Ep) 
for four different temperatures. At zero and low temperatures, 
in the magnetic phase, T < Tc, the Fermi level lies inside 
the hybridization gap for the up-band, but inside the conduc- 
tion band for down-spin, (k), which implies a semi-metal 
behavior in the magnetic phase. Then, an insulating phase is 
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FIG. 2: Plot of A|, A|, M and m as a function of temperature. Self- 
consistent solutions are obtained for the following parameters: Jx — 
0.8, Jh = -0.01 and Uc = 0.8 (in units of the half-bandwidth D). 



obtained for T > Tc when the Fermi level is inside the gap 
for both up and down spin directions and coincides with the 
energy of the f-level, Eq, (see. Fig. (3); and d]). Finally, when 
Act goes to zero, the hybridization gap closes and the system 
becomes metallic, but with a strong heavy fermion behavior 
due to the fact that the Fermi level coincides with the local- 
ized f-level, Eq. 

At last, from the quasiparticle spectrum @ one can es- 
timate the mass enhancement, which will be also spin- 
dependent lEHl : 



m 



1 



2a 



(Eo^ -Ep- A,)2 



(8) 



We can see in the Figl^that for both spin directions the mass 
enhancement increases with temperature and becomes dra- 
matic in the pure Kondo regime, when the magnetization goes 
to zero. As the denominator of (|8} goes to zero in the pure 
Kondo region the effective mass goes strictly to infinite. But 
in order to better visualize its behavior and particularly the 
second transition at Tk, we included a very small finite width 
in the /—level so obtaining extremely high values for the ef- 
fective mass in the pure Kondo region. Hence, two peaks are 
evident, one corresponds to the Curie temperature and the sec- 
ond, at higher temperatures, to the Kondo temperature. In the 
region of coexistence, the effective mass increases as a func- 
tion of temperature following the position of the Fermi level 
inside the E"_ (fc) band. Although simplified, this result pro- 
vides a qualitatively good explanation of the behavior of the 
electronic specific heat. 

In conclusion, we have formulated a mean field theory of 
coexistence of ferromagnetism and heavy fermion behavior 
in terms of self-consistent equations for the relevant fields. 
Within this treatment of the UKL model we obtain - for values 
of Jk between 0.7 and 1 .2 - a region where the order parame- 
ters. Act, the /—magnetization, M, and the c— magnetization, 
m, are different from zero, characterizing a coexistence be- 
tween the heavy fermion (Kondo) properties and ferromag- 
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FIG. 3: Temperature variation of tiie up- (solid line) and down- 
(dashed line) /-density of states, pfa{^~ Ef). Panel (a) corresponds 
to T = 0, (b) T = 0.02, (c) T = 0.06 and (d) T = 0.08. The pa- 
rameters are the same as in Fig|2| 



T = by varying Kondo coupling parameter Jk, but this 
time it is discontinuous and leads to a non-Kondo magneti- 
cally ordered state, as it can be observed in Fig^ We remark 
that this coexistence is mainly due to the presence of both one 
non-hybridized and one hybridized /-level (see FigO, and 
this ingredient is the key for the existence of the ferromagnet- 
ically ordered Kondo state. Thus, the proposed UKL model 
can explain the behavior of some uranium compounds where a 
ferromagnetically ordered state coexists with a heavy fermion 
Kondo state. The FM-Kondo coexistence arises from the fact 
that the localized spins are only partially screened, and conse- 
quently they keep their magnetic interaction at low tempera- 
tures. 
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MPIPKS, Dresden, and MPIDSO, Gottigen, for hospitality. 
The authors would like to acknowledge Brazilian agencies 
CNPq and FAPERGS for financial support. 




0,00 0,02 0,04 ^ 0,06 0,08 



FIG. 4: Temperature dependence of the effective mass enhancement. 
The parameter set is the same as in Fig|2| 

netic order. When increasing temperature the magnetic order 
disappears at the ordering temperature Tc while is still dif- 
ferent from zero, so indicating a Kondo or NFL regime. This 
is in agreement with experiments for uranium compounds 
where a Kondo behavior is observed also above the Curie tem- 
perature, like the examples we have quoted before. We remark 
also that the characteristic Kondo temperature, T^-, increases 
when the Kondo coupling value Jk also increases. Those 
results are completely different from the previous ones ob- 
tained for the normal KL model, where the Kondo effect and 
the magnetic order are in competition and cannot be present 
together The critical temperature Tc characterizes 

here a continuous quantum phase transition from the ferro- 
magnetic Kondo phase to the non-magnetic Kondo phase at 
higher temperatures. Another transition can be obtained at 
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